Controllable Andreev retroreflection and specular Andreev reflection in a 
four-terminal graphene-superconductor hybrid system 
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We report the investigation of electron transport through a four-terminal graphene- 
superconductor hybrid system. Due to the quantum interference of the reflected holes from two 
graphene-superconductor interfaces with phase difference 8, it is found that the specular Andreev 
reflection vanishes at 9 — while the Andreev retroreflection disappears at 6 = 71". This means that 
the retroreflection and specular reflection can be easily controlled and separated in this device. In 
addition, due to the diffraction effect in the narrow graphene nanoribbon, the reflected hole can exit 
from both graphene terminals. As the width of nanoribbon increases, the diffraction effect gradually 
disappears and the reflected hole eventually exits from a particular graphene terminal depending on 
the type of Andreev reflection. 
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Graphene, a single layer honeycomb lattice consist- 
ing of carbon atoms, has attracted considerable atten- 
tions in condensed matter community recently. [l|, [2], y, |4| 
The unique band structure of graphene with a linear 
dispersion relation near the Dirac-points leads to many 
peculiar properties, such as the low-energy Dirac-like 
quasi-particle dispersion relation and the relativistic-like 
behaviors. 0,13] Very recently, people begun to investi gate 
graphene-superconductor hybrid s vst ems . (El Irj M. M, lij. 1 1 0\ 
A unique and interesting phenomenon, the specular An- 
dreev reflection (different from the usual Andreev re- 
flection), was predicted to occur at the interface of the 
graphene and superconductor. [E| It was discovered fifty 
years ago.fllT] that near the interface of a conductor 
and superconductor an incident electron from the metal- 
lic side is retro-reflected as a hole and a Cooper pair 
is created in the superconductor, a process known as 
the Andreev reflection. When the bias is smaller than 
the superconductor gap, the conductance of the metal- 
superconductor hybrid device is mainly determined by 
the Andreev reflection. For the graphene-superconductor 
system, in addition to the Andreev retroreflection, an 
unusual Andreev reflection, the specular Andreev reflec- 
tion may occur, in which the direction of reflected hole 
is along the specular direction. [E| From the band struc- 
ture point of view, if electron-hole conversion is intra- 
band: both incident electron and reflected hole are from 
the same band (conduction or valence band), this corre- 
sponds to the usual Andreev retroreflection. The specu- 
lar Andreev reflection occurs if the electron-hole conver- 
sion is interband: the incident electron and reflected hole 
are, respectively, in the conduction and valence bands. 
Note that in two-terminal superconductor-graphene de- 
vice, both specular reflection and retroreflection occur. 
It is highly desirable to control and separate these An- 
dreev reflections experimentally. It is the purpose of this 



letter to achieve this goal. 

In this letter, we study a four-terminal graphene- 
superconductor device which consists of two supercon- 
ductor terminals with the phase difference and two 
graphene terminals (see Fig. la). By using the non- 
equilibrium Green function method, the current as well 
as the Andreev reflection coefficients are calculated. Our 
result shows that due to the quantum interference of re- 
flected holes from two superconductor terminals 2 and 
4 different Andreev reflection processes can be selected 
by tuning the phase difference 8. When = only the 
Andreev retroreflection occurs and the specular Andreev 
reflection is prohibited while for 6 = n only the spec- 
ular Andreev reflection occurs and retroreflection van- 
ishes. Therefore it is very easy to control the specular 
Andreev reflection and Andreev retroreflection by simply 
tuning the superconductor phase difference 9. In addi- 
tion, the direction of the reflected hole, which is along 
either the graphene terminal- 1 or terminal-3 depending 
on the type of Andreev reflection, can only be exhibited 
for large samples. When the sample size is comparable to 
the wavelength of reflected hole, however, the diffraction 
effect dominates so that the reflected hole can exit from 
both graphene terminals. 

The four-terminals device we considered consists of a 
zigzag edged graphene nanoribbon sandwiched by two 
superconductor terminals, as shown in the Fig. la. 
In the tight-binding representation, the Hamiltonian of 
the clean graphene nanoribbon is given bv[l3| Hg = 

T,ia E O a L a i<r + E<ii>a ta L a >> where 4r 0^) IS the 

creation (annihilation) operator at the site i. The 
on-site energy Eq is the reference energy for Dirac- 
point, which can be controlled experimentally by the 
gate voltage. Two superconductor terminals are repre- 
sented by BCS Hamiltonian, H Sa = Y, ka £kCl^ a C ka , a + 

E u (A a t7 ki , a C7_ kT , a + A* d kT ), where a - 2,4 



2 



FIG. 1: (color online) (a) is the schematic diagram for four- 
terminal graphene-superconductor device. In this diagram, 
the width of graphene nanoribbon is W = 6. (b) Tiu and 
TnA vs. the energy E for a three-terminal device with Eo = 
-0.5A and W = 25. 



is the index of the superconductor terminal and A a — 
Ae <" with the superconductor gap A and phase 9 a . 
The coupling between superconductor terminal a and 
graphene is described by Hxa = ta^ ia C a .,j(xi) + h.c. 
Here Xi is the horizonal position of the carbon atom i and 
Cka,cr-[lJ| So the total Hamilto- 



FIG. 2: Tua and T 13A vs. the energy E for W = 25. 



{{Hsa + Hto)- 



C a ^(x) = Y^k x ,ky u 

nian is H = H G + J2 a =2,4\ 

Using the Heisenberg equation of motion, [15| the cur- 
rent flowing from the graphene terminal 1 to the scatter- 
ing region is found to be 
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[(/i+ - h)Ti2 + (fi+ - U)T U 



2e f dE 

(fl + - f3-)Ti3A + (/l+ - fl-)TuA 
(fl+ - /3+)Tl3]- 



(1) 



where f a ±(E) = l/{exp[(E T eV a )/k B T] + 1} and 
f 2 (E) = U(E) = l/{exp(£;/fc B T) + l} are the Fermi dis- 
tribution with the bias V a . Here we set the bias of two su- 
perconductor terminals be zero (V2 = V4 = 0). In Eq.(l), 
T 13 (E) = Tr{r 1TT G!f T r3 TT G° T } and T 12[14) (E) = 
Tr{ri||[G r r 2 (4)G a ]||} are the normal transmission co- 
efficients from the terminal 1 to the terminal 3, 2, and 
4, respectively. T nA (E) = Tr{ri TT G^r lu GJ T } and 
Ti3A(E) = Tr{r 1TT G^r 3U G| T } are the Andreev reflec- 
tion coefficients for the incident electron coming from the 
terminal 1 with the hole Andreev reflected to the termi- 
nal 1 (Tua) or terminal 3 (T 13j 4). Here the subscripts 
T1\ t-lo IT; and II represent the 11, 12, 21, and 22 ma- 
trix elements in Nambu subspace. The linewidth function 
T a (E) is defined as T a (E) = i[W a - (W a )^] and G r[a) (E) 
are the retarded (advanced) Green functions of central 
region m Nambu representation. G r (E) = G a ^(E) = 
(EI — H c — J2 a =i 234 ^q) _1 with the Hamiltonian H c 
of the central region labeled by a rectangular area in 
Fig. 1(a). S^(-E) is the retarded self-energy due to the 
coupling to the terminal a. Y, r atj (E) = tg^ ij (E)t, 
where g r a ^ (E) is the surface Green function of termi- 
nal a. For the graphene terminal 1 and 3, we have to 
numerically calculate their surface Green function.flfl] 
while for superconductor terminal 2 and 4, the sur- 
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face Green function 

A*~/E \ >[3| where p is the normal den- 

sity of states, Jo[kp(xi — Xj)} is the 0-th order Besscl 
function with the Fermi wave vector fcp, and (5(E) — 
—iE/\/A 2 - E 2 for \E\ < A and (3(E) = \E\/\/E 2 — A 2 
for \E\ > A. In numerical calculations, we set the hop- 
ping energy t = 2.75eV and the length of C-C bond 
ao = 0.142nm as in a real graphene sample. The su- 
perconductor gap A is set to be A = lmeV and the 
Fermi wave- vector hp = lA 

We first study a three-terminal device by decoupling 
one of the superconductor terminal (2 or 4). Fig. lb shows 
the Andreev reflection coefficients Tua and Ti^a as a 
function of incident electron energies E. It can be seen 
that Tn a and Ti^a are quite large when the energy E is 
within the gap (\E\ < |A|) and exhibit peaks at the Dirac 
points E = ±Eq and the gap edge E = ±A. 17[ Similar 
to the usual normal-superconductor junction Tua and 
Ti 3 a decay quickly when E is outside of the gap.fljj. 
Note that when \E\ < \E \, the incident electron and re- 
flected hole are in the same band (see Fig.3f) leading to 
the usual Andreev retroreflection. On the other hand, 
for \E\ > \Eq\, the incident electron and reflected hole 
are, respectively, in the conduction and valence bands 
(see Fig.3f) giving rise to the specular Andreev reflec- 
tion. The above results show that both retroreflection 
and specular reflection occur with large amplitudes in 
the three-terminal device (with only one superconductor 
terminal) . 

Next, we focus on the four-terminal device. Fig. 2 
shows Tn a and Ti^a versus the energy E for two dif- 
ferent superconductor phase differences 9 = 9 2 — 0± = 
and 7T. When 9 — 0, Tua and Ti 3 a are zero for \E\ > \E a \ 
but quite large for \E\ < \E \ (see Fig.2a and 2b). This 
means that at 9 = only the retroreflection occurs and 
the specular reflection is prohibited. On the other hand, 
when 9 — 7r, the situation reverses: Tua and Ti%a are 



3 



FIG. 3: (color online) (a)-(d) is T 13A vs. the width W and (e) 
is Tiia vs. W with the parameters Eo = — 0.7A and 9 — 0. 
(f) is the schematic view of the Andreev retroreflection and 
specular Andreev reflection. 



zero for \E\ < \E \ and quite large when \E\ > \E \ (see 
Fig. 2c and 2d). Hence when 9 — ir, the retroreflection is 
prohibited and only specular reflection occurs. Experi- 
mentally, the phase difference 9 can be tuned by varying 
the super-current between two superconductor terminals. 
So the present four-terminal device gives us a handle to 
experimentally control and select the Andreev retrore- 
flection and specular Andreev reflection. 

Now we explain why the retroreflection disappears at 
the phase difference 9 = tt while the specular reflection 
vanishes at 9 — 0. In the four-terminal device with two 
graphene-superconductor interfaces, two Andreev reflec- 
tions from each interface contribute coherently to the re- 
sultant Andreev reflection coefficient. Depending on the 
phase carried by each Andreev reflection, the interference 
can either be constructive or destructive. For the retrore- 
flection, each reflected hole carries a phase factor, j 1 81] 9 a 
of the corresponding superconductor terminal, leading 
to a total Andreev reflection coefficient proportional to 
\e 102 + e l9i \ 2 = 1 1 + e te \ 2 , whose value reaches the maxi- 
mum at 9 = and minimum at 9 — tt. So the Andreev 
retroreflection disappear at 9 = tt due to the destructive 
interference. However, for the specular Andreev reflec- 
tion, in addition to the phase difference 9, an extra phase 
7r is acquired due to the reflection between two interfaces 
when the incident electron and reflected hole involves dif- 
ferent energy bands. The origin of this extra phase tt is 
similar to the tt junction of the superconductor-graphene- 
superconductor device, [6[ where a super-current of form 
I = I c sin(# + 7t) was found. Due to this extra phase 
the total Andreev reflection coefficient is proportion to 
|1 + e ! ( e ' +7r )| 2 , whose value is zero at 9 = resulting a 
vanishing specular Andreev reflection at 9 = 0. 

Since the reflected hole from the Andreev retroreflec- 



FIG. 4: (color online) T 11A and T 13A vs. W for E = -0.7A 
and 8 — tt. 



tion (the specular Andreev reflection) is along the retro- 
reflected (the specular reflected) direction as shown in 
Fig. la, Andreev reflection coefficient T 13 a (Tua) should 
be zero at 9 — (tt). This does not agree with what we 
have obtained in Fig. 2. We attribute this phenomenon 
to the diffraction effect of the reflected hole in the small 
device. To verify this statement, we have studied the size 
dependence of the Andreev reflection coefficient. Fig. 3a- 
3e show Ti3A and Tua versus the width W for 9 = 0, 
in which only the retroreflection occurs. With the in- 
crease of the width W, Ti 3 a oscillates and decays to 
zero while Tua increases and saturates at large W. This 
clearly indicates that the reflected hole exits only from 
the terminal-1 at large width W. In addition, from the 
standing wave in the semi-infinite graphene ribbon, the 
wavelength of the reflected hole has been calculated. We 
found that when the graphene-ribbon width W is in sev- 
eral tens this wavelength is on the same order of the de- 
vice size. So for this width W (several tens), the diffrac- 
tion effect is significant. As a result, the reflected hole 
can exit from both terminal 1 and 3 leading to non-zero 
values of both Tua and Ti^a- 

We now examine the size dependence of the specular 
Andreev reflection when 9 = tt. We expect that when the 
width W of the graphene ribbon increases the reflected 
hole should go to the terminal 3 due to the fact that 
the diffraction effect disappears at large W . Indeed, our 
numerical result confirms this. From Fig. 4 we see that 
Tua decays to zero and Ti 3 a saturates at large W. 

From the above discussion, we see that when the wave- 
length of the reflected hole is comparable to the width 
of graphene nanoribbon, its direction can not be used 
to distinguish the Andreev retroreflection and specular 
Andreev reflection. Nevertheless, these two kinds of An- 
dreev reflection can manifest their difference in the four- 
terminal device by tuning the superconductor phase as 
demonstrated above. This also shows that the electron- 
hole conversion mechanism, i.e., interband or intraband 
conversion, is the fundamental origin for these two kinds 
of Andreev reflections. 

Up to now, we only considered the clean graphene rib- 
bon at 9 = or tt. In the presence of the weak impu- 
rity disorder, [19( all the results still remain, except that 
the boundary of two kinds of Andreev reflection slightly 
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FIG. 5: Upper panel: Tha (a) and T 13 a (b) vs. E for different 
8. Lower panel: the conductance G of terminal 1 as a function 
of bias V. The parameters are Eo = — 0.5A and W = 50. 



smeared. In addition, for other 9, both types of Andreev 
reflection may occur (see Fig. 5a and 5b) due to the in- 
complete destructive interference. With the variation of 
9 from to tt, the specular reflection (Tha and T\za a t 
\E\ > \Eq\) gradually increases from zero to the maxi- 
mum value while the retroreflection (Tha and T\sa at 
\E\ < \Eq\) gradually decreases from the maximum value 
to zero. 

Finally, we investigate the differential conductance 
G = dh/dV at zero temperature and discuss the ex- 
perimental feasibility. By setting the bias V\ = V3 = V 
and V2 = Va = 0, the direct tunneling T13 from the 
terminal 1 to the terminal 3 does not contribute to 
the current. The differential conductance is given by 
G{V) = {2Tua + 2T 13A + T12 + T u )2e 2 /h. At small 
bias \eV\ < A, the normal tunneling processes from the 
terminal 1 to two superconductor terminal 2, 4 (T 12 and 
T 14 ) are forbidden, so the conductance G{V) is directly 
related to the Andreev reflection coefficient T\\a +T%sa- 
Fig. 5c shows the conductance G versus the bias V at 
different 9. For 9 — 0, the specular Andreev reflection 
vanishes, so G is zero at \Eq\ < \eV\ < A. On the other 
hand, for 9 = tt, the Andreev retroreflection disappears 
leading to G — at \eV\ < \Eq\. Upon varying 9 from 
to 7r, the conductance G at \eV\ < \E \ drops to zero 
while G at \Eq\ < \eV\ < A gradually increases from 
0. For the bias \eV\ > A, G is always large and weakly 
depends on the phase difference 9 because of the contri- 
bution of the normal tunneling process. Note that exper- 
imentally the graphene nanoribbonsj^] have been real- 
ized with the sup erconductor leads attached to graphene 
nanoribbons.[10| So the proposed device is within the 
reach of the present technology and is feasible experi- 
mentally. 

In conclusion, the interplay of two kinds of Andreev 
reflections in a four-terminal graphene-superconductor 
hybrid device was investigated. It was found that the 



Andreev retroreflection and specular Andreev reflection 
can be tuned in this system due to quantum interfer- 
ence. When the superconductor phase difference 9 = 0, 
the specular Andreev reflection is prohibited and only 
the Andreev retroreflection occurs. However, for 9 = tt, 
the Andreev retroreflection is suppressed and only the 
specular Andreev reflection occurs. In addition, in the 
narrow graphene nanoribbon with its size comparable to 
the wavelength of the reflected hole, the diffraction ef- 
fect occurs. Then the reflected hole can exit from both 
graphene terminals. On the other hand, for large sam- 
ples, the diffraction effect disappears and the reflected 
hole can only traverse to a particular terminal depend- 
ing on the kind of Andreev reflections. 
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